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Abstract. The phase-integral approximation devised by Froman and Froman, is used for 
computing cosmological perturbations in the quartic chaotic inflationary model. The phase- 
integral formulas for the scalar power spectrum are explicitly obtained up to fifth order of the 
phase-integral approximation. As in previous reports [1—3], we point out that the accuracy 
of the phase-integral approximation compares favorably with the numerical results and those 
obtained using the slow-roll and uniform approximation methods. 
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1 Introduction 

Chaotic inflationary models like V oc (fi p have been studied and analized with the WMAP 
data in recent years. The recent observations reported by WMAP favors inflation, and the 
quartic chaotic model, with p = 4 [4, 5] still remains ruled out by this data [6, 7]. Ramirez 
et al. [8] point out that if the scalar power spectrum is considered not scale invariant, then 
quartic models would not be excluded by the recent cosmological data. Ramirez et al. [8] 
consider a new scenario where the observable modes cross the horizon near the beginning of 
inflation. 

In order to compare with observations, we should be able to obtain very accurate results 
for the predicted power spectrum of primordial perturbations for a variety of inflationary sce- 
narios. In general, most of the inflationary models are not exactly solvable and approximate 
or numerical methods are mandatory in the computation of the scalar and power spectra. 
Traditionally, the method of approximation applied in inflationary cosmology is the slow- 
roll approximation [9], which produces reliable results in inflationary models with smooth 
potentials, but cannot be improved on a simple way beyond the leading order. Recently, 
some authors have applied alternative approximations, such as the WKB method with the 
Langer modification [10-12], the Green function method [13], and the improved WKB method 



Habib et al [16-18] have successfully applied the uniform approximation method in the 
calculation of the scalar and tensor power spectra and the corresponding spectral indices 
for the quadratic and quartic chaotic inflationary models, showing that the the uniform 
approximation gives more accurate results than the slow-roll approximation. Casadio et al 
[19] have applied the method of comparison equation to study cosmological perturbations 
during inflation. The comparison method is based on the uniform approximation proposed 
by Dingle [20] and Miller [21] and thoroughly discussed by Berry and Mount [22]. 

Recently [1-3], the phase integral method has been used in the calculation of the scalar 
and tensor power spectra and the spectral indexes for the power-law and quadratic chaotic 
model, obtaining better results than those derived using slow-roll [9], WKB and the improved 
uniform approximation methods [16-18]. 



[14, 15]. 
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Since the 7- year WMAP data does not show evidence for tensor modes [23, 24], in 
this article we use the phase-integral approximation up to fifth order to calculate the scalar 
power spectrum for the quartic chaotic inflationary model. We have calculated the tensor 
power spectrum and since the results are very similar to those obtained in the scalar case, 
we decided to skip their presentation in this article. f 

The article is structured as follows: In Sec. 2 we numerically solve the equation gov- 
erning the scalar perturbations. In Sec. 3 we apply the phase-integral approximation to the 
quartic chaotic inflationary model. We present the improved uniform approximation method, 
the slow-roll approximation and the numerical solution to the perturbation equation. The 
comparison among different methods is presented in Sec. 4. Finally in Sec. 5 we present our 
conclusions. 



2 Equation for the perturbations 

In an inflationary universe driven by a scalar field, the equations of motion for the inflaton 
4> and the Hubble parameter H are given by [25] 



+ 3H(j) 



dV(<f>) 



(2.1) 



H 



1 



V(<f>) + 



(2.2) 



with V(4>) = \\<f) A and the constant A, can be fixed with the help of the amplitude of the 
density fluctuations observed by WMAP. A fitting to the observational data requires that 
A ~ 1CP 13 [18]. Using Eq. (2.1) and Eq. (2.2) we obtain that, in the slow-roll approximation, 
the expansion factor a(t) and the inflaton field (p(t) are 



a ~ a,- 



1 — exp 




Mpit 



4>{t) ~ ^exp [ -2\l -Mpii I , 



(2.3) 
(2.4) 



where a.; and <f>i are integration constants corresponding to the initial values of the scale 
factor and the inflaton. 

The application of the the phase-integral approximation for solving the equations gov- 
erning the scalar perturbations requires the knowledge of the evolution of the expansion 
factor a and the scalar field (j), and the computation of analytic expression of them would 
highly simplify the calculations. In order to find (j) and a as a function of time, we solve the 
coupled system of differential equations (2.1)-(2.2) with the sixth-order Runge-Kutta method 
[26] and fit the numerical data with the software gnuplot following the form of the slow-roll 
solutions (2.3,2.4). The integration is performed in the physical time t using the sixth-order 
Runge-Kutta method [26]. The initial value for the inflaton is 4>i = 21.9 Mpi. The initial 
condition for the speed of the scalar field <j>i is obtained from the slow-roll solution (2.4). The 
initial value for aj is chosen as unity and A = 1.75 x 1CP 13 . These initial conditions guarantee 



sufficient inflation. We find that the inflation ends at tf 
corresponding to 59.98 e-folds before the scalar field starts to oscillate. 



4.62 x 10 6 M" 1 



1.25 x 10 



-36 c 
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The equation governing the evolution of the scalar perturbations in the quartic chaotic 
inflationary model is given by 



a 1 

Uk + -Uk + -~ 
a a z 



k 



(azs + az s ) a 



Uk = 0. 



(2.5) 



In order to apply the phase-integral approximation we need the expressions for a(t) and 
4>{t) which have been derived fitting the numerically solution of Eqs. (2.1) and (2.2). We 
have obtained the following analytic expressions 



a(t) = di exp {a [l — exp ( — /3t)j } , 
4>(t) = 0j exp (-7*) , 



(2.6) 
(2.7) 



where a = 60.2839, (3 = 9.6370 x 10~ 7 M P1 , and 7 = 4.8305 x 10~ 7 M P i and, the function 
zs(t) has he form 



zs(t) = -ai4>i^= exp {a [l - exp (-/3t)] + fit - jt\ . 



(2i 



Since we want to eliminate the term iik in Eq. (2.5), we introduce the change of variable 
u kif) = Uk r- > i obtaining that Uk satisfies the differential equation: 



U k + Rs(k,t)U k = 0, 



(2.9) 



with 



Rs(k,t) = 4t 



(azs + az's) a 



2S 



+ 4^2 (« 2 " 2a «) 



(2.10) 



where U(k) satisfies the asymptotic conditions 



Uk ->■ A k \/a(t)zs(t), kt^t 00, 
U k - 



a(t) 



2k 



exp [—ikrj(t)], kt — > 0. 



(2.11) 
(2.12) 



We now proceed to write the explicit equations for quartic chaotic inflation. From Eq. 
(2.10), with Eq. (2.6) and Eq. (2.7) we obtain 



k 2 9 
Rs(k,t) = — exp{-2a[l-exp(-/3t)]}--a 2 /3 2 exp(-2/3t) (2.13) 

a i 4 
- ^(^-27)exp(-^)-^- 7 ) 2 . 

In order to apply the asymptotic condition (2.12), we use the relationship between rj 
and the physical time t, which is given in terms of the integral exponential function Ei(z) 
[27] as 
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rj(t) = — = exp (—a) {Ei [q exp (— /3to)] — Ei [a exp (— f3t)] } , 



(2.14) 



where t{ = 10 7 Mpj 1 , so that r\ is zero at the end of the inflationary epoch. The dependence 
of r\ on t is shown in Fig. 1. We can observe that, as 



- krj 
-kr) - 



■> 

oo 



kt 
kt 



oo, 
0. 



(2.15) 
(2.16) 
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Figure 1. Behavior of r\ as a function of f for the chaotic inflationary \ Ac/) 4 model. 



Eq. (2.9), where Rs(k,t) is given by Eq. (2.13) , does not possess exact analytic 
solution. In order to solve the differential equations governing the scalar perturbations in the 
physical time t, we use the fifth-order phase integral approximation [28] and compare this 
results with those obtained using the slow-roll and the improved uniform approximation. 

Once the mode equation for scalar perturbations is solved for different momenta k, the 
power spectra for scalar modes is given by the expression 



P s (k) = lim 

k v ^o 2ir 2 



Uk(rj) 



z{rj) 



(2.17) 
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The spectral indice is defined as 



nW = 1+ *>EMm. (2.i8) 

3 Methodology 

3.1 Phase-integral approximation 

In order to solve Eq. (2.9) with the help of the phase-integral approximation [28], we choose 
the following base functions Q$ for the scalar perturbations 



Ql(M) = Rs{k,t), (3.1) 

where R§(k,t) is given by Eq. (2.13). Using this selection, the phase-integral approximation 
is valid as kt — > oo, limit where we should impose the condition (2.11), where the validity 
condition fi <C 1 holds. The selection, given in Eq. (3.1), makes the first order phase-integral 
approximation coincides with the WKB solution. The bases functions Qs(k, t) possess turning 
points t rc t = v s = 117726 Mpj 1 for the mode k = 0.05 Mpc -1 . The turning point represents 
the horizon. There are two ranges where to define the solution. To the left of the turning 
point < t < t rct we have the classically permitted region Qg(k,t) > and to the right of 
the turning point t > t Tet corresponding to the classically forbidden region Qg(k,t) < 0, such 
as it is shown in Figs 2(a). 

The mode k equations for the scalar perturbations (2.9) in the phase-integral approxi- 
mation has two solutions: For < t < t TCt 



«?(t) 



+ 



C2 



q- 1/2 (k,t) 


cos 


\ws{k,t)\ - j 


?s 1/a (M) 


COS 


|w S (M)| + | 



(3.2) 



For t > t rct 



«?(*) 



ci 



+ 



C2 



q- 1/2 (k,t) 



exp 



\us(k,t)\] 



(3.3) 



1/2 (k,t) exp[|w s (*!,t)|] 



The phase- integral approximation solutions are of the order 2N + 1, so N = 2 means 
that we have the approximation up to fifth order and qs(k,t) can be expanded in the form 



2 

«s(M) = £WM)Qs(M) = [Y Qs (k,t) +Y 2s (k,t) +Y 4s (k,t)]Q s (k,t). (3.4) 

n=0 

In order to compute qs(k,t), we compute Y2 s (k,t), Y$ s (k,t), and the required functions 
eo s (k,t), e2 s (k,t). The expression (3.4) gives a fifth-order approximation for qs(k,t). In 
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Figure 2. (a) Behavior of Qa(k, t). (b) Contour of integration T vs (t) for < t < ug- (c) Contour 
of integration T vs (t) for t > i/g- The contour of integration has been done in the complex z plane. 
For this function we have a Riemann surface with two sheets where the function is single- valued. In 
the contour of integration the solid line indicates the first Riemman sheet whereas that the second 
Ricmman sheeet is indicated for the dashed line [29]. 
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order to compute us(k,t) we make a contour integration following the path indicated in Fig. 
2(b)-(c). 



LO S (k,t) = U 0s (k,t) + ^2u) 2 n s (k,t), 
n=l 

i 2 r 

Qs(k,t)dt + -J2 / Y 2ns (k,t)Q s (k,t)dt, 

n=l ^ r ^s 
2 

/ Q s (k,t)dt + ^^2 I f2n s (k,t)dt, 



/ 



where 



f 2ns (k,t) = Y 2ns (k,t)Q s (k,t). 
The functions f 2ns {k,t) have the following functional dependence: 



(3.5) 
(3.6) 
(3.7) 

(3.8) 



/2s (M) 

h s (k,t) 



A(k,t)(t-v s ) 
B(k,t)(t-v s ) 



-5/2 
-11/2 



(3.9) 
(3.10) 
(3.11) 



where A(k, t) and B(k, t) are regular at vs ■ With the help of the functions (3.9) we compute 
the integrals for uj 2n up to N = 4 using the contour indicated in Fig. 2(b)-(c). The expressions 
for U2n permit one to obtain the fifth-order phase integral approximation of the solution to 
the equation for scalar perturbations (2.9). The constants c\ and c 2 are obtained using the 
limit kt — > of the solutions on the left side of the turning point (3.2), and are given by the 
expressions 



ci = -ic 2 , 



c 2 



V2 



-i fcr?(0)+L»o s 0,0) 



(3.12) 
(3.13) 



In order to compute the scalar power spectrum, we need to calculate the limit as kt — > oo 
of the growing part of the solutions on the right side of the turning point given by Eq. (3.3) 
for scalar perturbations. 



«?(*) 



*s(*) 



(3.14) 



Fs(k)= lim 
3.2 Improved Uniform approximation 

We want to obtain an approximate solution to the differential equation (2.9) in the range 
where Qg(k,t) have a simple root at t ret = vs, so that Qg(k,t) > for < t < t TCt and 
Qs(k,t) < for t > t TCt as depicted in Fig. 2(a). Using the uniform approximation method 
[1-3, 16, 30], we obtain that for < t < t TCt we have 
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U k (k,t) 
[pi(k,t)f 2 



Qi(k,t) 



1/4 



{CiAi[-pi(k,t)] + C 2 Bi[-pi(k,t)}} 



[Ql(k,t)] 1/2 at 



(3.15) 



(3.16) 



where C\ and C 2 are two constants to be determined with the help of the boundary conditions 
(2.12). For t > t rct 



U k (k,t) 



-Pr(k,t) 

Ql(k,t) 



i/i 



{C 1 A i [p T (k,t)]+C 2 B i [p r (k,t)}} 



[Pr(k,t)f 2 = f [-Qi( k ,t)] 1/2 at, 

J trrt 



For the computation of the power spectrum we need to take the limit kt 
solution (3.17). In this limit we have 



(3.17) 

(3.18) 
oo of the 



£= [-Ql(k, t)] " 1/2 { \ exp (- £ [-Ql(k, t)] 1/2 at 



+ iexp(£ [-Q 2 s (k,t)} 1/2 dt^y 



(3.19) 



where C is a phase factor. Notice that Eq. (3.19) is identical to Eq. (3.3) obtained in 
the first-order phase-integral approximation. Using Eqs. (2.6), (2.8) and the growing part 
of the solutions (3.19), one can compute the scalar and power spectrum using the uniform 
approximation method, 



bV ; -kt^oo2lT 2 Z S (t) 



(3.20) 



We using the second-order improved uniform approximation for the power spectrum [18], 

(3.21) 



P s (k) = P s (k) [T*(v s )}, 
where v§ is the turning point for the scalar power spectrum and 

139 



r»EEi + -L + 



+ 



(3.22) 



I2v 288v 2 51840f 3 
3.3 Slow-roll approximation 

The scalar power spectra in the slow-roll approximation to second-order is given by the 
expression [13, 31] 



P^(k) 



7tT< 



1 + (46 - 2)ei + 2b5i + ( 4b 2 - 23 + — J e( + ( 3b 2 + 2b - 22 + ) €l 5i 



29tt' 



5vr 2 



+ [3^-4+— St + l-b< + -)6 2 



it 



( H_\ 2 (H^ 



I 2?r J 



(3.23) 



k=aH 
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where b ~ 0.7296 is the Euler constant. The spectral index in the slow-roll approximation to 
second-order is 

n| r (fc) ~ 1 - 4ei - 25i + (8c - 8)e? + (10c - 6)e 1 5 1 . (3.24) 

The parameters used in the equations (3.23) and (3.24), express in terms of the potential 
V(4>), are given by [13] 

(f ) 2 (» + > ~ > ~ & + " I* - 1 14 ) ■ < 326) 

fi - 5^1 - |k? + \um + \ul - jjt^y, + ^t/.v, 2 + ii/^j, (3.27) 

1 2 4 1 1 3 23 

*i ^ 2^ " Fl " 3^ + 3 UlVl ~ 3 V ? ~ 3 V2 + 2 U * ~ + ~9 UlV ? (3-28) 

h * Ul - |t7iVi + If + V 2 - + ^C/fVi - ^E/il? + (3.29) 
2 bo 2 3 

" 1^2 + 2^2 + ^3. 

With, 



I7i = 




(3.30) 


Vi = 


V" 
V ' 


(3.31) 


V 2 = 


y'V" 


(3.32) 


v 3 = 


(V'f V"" 
V 3 ' 


(3.33) 



The expressions (3.25)-(3.33) depend explicitly on <fi. In order to compute the scalar 
power spectrum we need to obtain the dependence on the variable k. For a given value of k 
(0.0001 Mpc^ 1 < k < 15 Mpc" 1 ) we obtain t from the relation k = aH. Thus, for each k one 
obtains a value of t that we substitute into Eq. (2.4). 

3.4 Numerical solution 

We integrate on the physical time t the equations (2.9) governing the scalar perturbations 
using the predictor-corrector Adams method of order 12 [26], and solve two differential equa- 
tions, one for the real part and another for the imaginary part Uk- Two initial conditions are 
needed in each case Uk(U), U k (ti), which can be obtained from the third-order phase-integral 
approximation. We start the numerical integration at ti calculated at 25 oscillations before 
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reaching the turning point t rct [32]. We call this procedure ICs phi3. Fig. 3 compare the 
numerical solution with the fifth-order phase-integral approximation for Re(«fc). Figures are 
plotted against the number of e-folds N. The solid line corresponds to the numerical solution 
(ICs phi3), the dashed line corresponds to the fifth-order phase- integral approximation. The 
turning point vs is indicated with an arrow. We stop the numerical computation of Ps(k) at 
t = 5.00 x 10 5 Mp 1 1 , after the mode leaves the horizon, where Uk/zs is approximately con- 
stant. Notice that the expressions for fitting (2.6) and (2.7) are valid in the aforementioned 
time scales, therefore we can use them for computing the scalar Ps(k) power spectrum. 
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Figure 3. Re(ufc) versus the number of e-folds for the chaotic inflationary j\(j) A model. Solid line: 
numerical result (ICs phi3); dashed line: fifth-order phase- integral approximation. 



4 Results 



For the chaotic \\<i) i inflationary model, we want to compare the scalar power spectra for 
different values of k calculated using the third and fifth-order phase-integral approximation 
with the numerical result (ICs phi3), the first and second-order slow- roll approximation and 
the first and second-order uniform approximation method. First we analyze the results for 
the scalar Ps(k) power spectrum shown in Fig. 4. 

Table 1 shows the value of Ps(k) and the relative error Ps(k) using each method of 
approximation at the WMAP pivot scale [18]. It can be observed that the best value is 
obtained with the fifth-order phase-integral approximation. It should be noticed that the 
slow-roll approximation works well since the parameters e\, ei and 5 n are small. 
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Figure 4. Ps(k) for the chaotic inflationary \\<j) A model. Thin solid line: numerical result (ICs phi3); 
dot-dashed line: third-order phase-integral approximation; dashed line: fifth-order phase-integral 
approximation; thick solid line: first-order phase-integral approximation, WKB and first-order uni- 
form approximation; dashed double-dots line: second-order improved uniform approximation, double- 
dashed dot line: second-order slow-roll approximation; dotted line: first-order slow-roll approximation. 
The inset is an enlargement of the figure. 



Method 


P s (k) x 1(T 9 


rel. error Ps(k)\ 


Numerical 


1.7995 




Third-order phase-integral approximation 


1.8021 


0.1407% 


Fifth-order phase-integral approximation 


1.7993 


0.0149% 


Leading order of slow-roll approximation 


1.7889 


0.5866% 


First-order slow-roll approximation 


1.7952 


0.2394% 


Second-order slow-roll approximation 


1.7953 


0.2324% 


First-order phase integral approximation, WKB, 


1.6181 


10.0851 % 


and first-order uniform approximation 






Second-order improved uniform approximation 


1.7857 


0.7694 % 



Table 1. Value of Ps(k) and its relative error obtained with different approximation methods for the 
chaotic inflationary model j\<f) 4 for the mode k = 0.05 Mpc -1 . 

The first-order phase integral approximation, the WKB and the first-order uniform 
approximation give the same result, and deviate from the numerical result in 10.0851%. 
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The second-order-improved uniform approximation gives an error of 0.7694%. With the 
leading, first and second-order slow-roll approximation we have an error for Ps(k) of 3> 100 %, 
0.2394% and 0.2324%, respectively. Using the third-order phase-integral approximation the 
error gives 0.1407%, whereas the fifth-order phase-integral, due to the convergence of the 
method, reduces to 0.0149% for Ps(k). Note that rel. error Ps(k) with srl ~ rcl. error Ps(k) 
with sr2, which is unexpected according to the results obtained previously for the power-law 
inflationary model [17]. The computation time is different for each method, for the scalar 
power spectrum we have the following times: numerical ~ several hours, pi3 « few minutes, 
pi5 ~ several minutes, srO ~ few seconds, srl ~ few seconds, sr2 « few seconds, WKB, ual, 
pil « few seconds, ua2 ps few seconds. Fig. 5 shows the results for the spectral indice ns(k). 




935 1 1 i i i 

' 0.0001 0.001 0.01 0.1 1 10 

k[Mpc 1 ] 

Figure 5. ns(fc) for the chaotic inflationary \X(j) 4 model. Thin solid line: numerical result (ICs phi3) 
and fifth-order phase-integral approximation; dot-dashed line: third-order phase-integral approxima- 
tion; thick solid line: first-order phase-integral approximation, WKB and first-order uniform approx- 
imation; dashed double-dots line: second-order improved uniform approximation; double-dashed dot 
line: second-order slow-roll approximation; dotted line: first-order slow-roll approximation. The inset 
is an enlargement of the figure. 

Table 2 shows the value of ns(k) and the relative error ns(k) using each method of 
approximation at the WMAP pivot scale [18]. The first-order phase integral approximation, 
the WKB and the first-order uniform approximation give the same result, and deviate from 
the numerical result in 0.007239 %. The second-order-improved uniform approximation gives 
an error of 0.003909%. With the leading, first and second-order slow-roll approximation we 
have an error for n s (k) of 6.0692239 %, 0.116637% and 0.0012000%, respectively. Using the 
third-order phase-integral approximation the error gives 0.000157%, whereas the fifth-order 
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phase-integral, due to the convergence of the method, reduces to 0.000026% for ns(k). Note 
that the first and second-order slow-roll approximation gives a better results for the spectral 
indice ns(k) than the power spectrum Ps(k). 



Method 


n s (k) 


rel. error ns(fc) 


Numerical 


0.942780 






Third-order phase-integral approximation 


0.942779 


0.000157 


7o 


Fifth-order phase-integral approximation 


0.942781 


0.000026 


7o 


Leading order of slow-roll approximation 


1 


6.069239 


7o 


First-order slow-roll approximation 


0.943880 


0.116637 


7o 


Second-order slow-roll approximation 


0.942894 


0.012000 


7o 


First-order phase integral approximation, WKB, 


0.942849 


0.007239 


7o 


and first-order uniform approximation 








Second-order improved uniform approximation 


0.942817 


0.003909 


7o 



Table 2. Value of ns(fc) and its relative error obtained with different approximation methods for the 
chaotic inflationary model jXtp 4 for the mode k = 0.05 Mpc -1 . 

The relative errors with respect to the numerical result have been obtained using the 
expressions: 



rel. error P s (k) = — x 100, 



pnum 



(*) 



rel. error n s (k) = ^ ^J m(fc) S > ] x 100. 



(4.1) 
(4.2) 



5 Conclusion 

As in our previous work [1-3], we show that the phase-integral approximation is an efficient 
method to calculate the scalar (or tensor) power spectrum and the spectral index in several 
models of inflation. Our results indicate that the scalar power spectrum is almost scale invari- 
ant since we have considered the usual scenario where the observable modes cross the horizon 
before the inflation ends, unlike the proposal presented by Ramirez et. al [8]. Therefore, the 
good agreement between the numerical results and those obtained with the phase-integral 
approximation shows that the phase integral method is a very useful approximation tool for 
computing the scalar and tensor the power spectra in a wide range of inflationar scenarios. 
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